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in  real  materials  based  on  the  identification  of  the  d-valence  as  a  critical  parameter  controlling  the  transition,  and  identifying  that  this  finding 
is  insensitive  to  the  choice  of  orbital  provided  it  is  well  localized.  Secondly,  we  discovered  a  new  mechanism  for  correlation-driven  metal- 
insulator  transitions  (the  site-selective  Mott  transition)  and  demonstrated  that  this  mechanism  accounts  for  the  main  features  of  the 
previously  mysterious  metal-insulator  phase  diagram  of  the  rare-earth  nickelates.  Thirdly,  we  demonstrated  the  ability  to  design  new  Mott 
insulators  via  heterostructuring.  Finally,  we  developed  a  fully  charge  self-consistent  implementation  of  DFT-l-DMFT,  which  allows  us  to 
compute  total  energies.  This  major  technical  development  allowed  us  to  compute  the  pressure  vs.  rare  earth  phase  diagram  in  the  nickelates. 
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This  is  the  final  progress  report  for  Electronic  Structure  at  Oxide  Interfaces,  contract 
number  W911NF-09- 1-0345  and  proposal  number  ARO  56032PH,  covering  the  period  Septem¬ 
ber  1st  2011  to  Aug  31,  2013.  The  broad  goal  of  the  project  is  the  identification  of  the  factors 
controlling  the  nontrivial  many-body  physics  in  artificially  nanostructured  materials  based 
on  transition  metal  oxides,  thereby  enabling  the  design  of  materials  with  desired  correlated 
electron  properties  such  as  ferromagnetism  with  a  high  Curie  temperature,  high  transition 
temperature  superconductivity  and  metal  insulator  transitions.  The  means  of  achieving  this 
goal  is  the  implementation  and  development  of  the  hybrid  methodology  of  density  functional 
theory  and  dynamical  mean-field  theory  (DFT-I-DMFT).  A  number  of  important  milestones 
were  achieved,  among  which  are: 

•  a  vastly  improved  understanding  of  the  physics  of  the  correlation-driven  metal-insulator 
transition  in  real  materials  based  on  the  identification  of  the  d-valence  as  a  critical 
parameter  controlling  the  transition,  and  identifying  that  this  finding  is  insensitive  to 
the  choice  of  orbital  provided  it  is  well  localized. 

•  the  discovery  of  a  new  mechanism  for  correlation-driven  metal-insulator  transitions 
(the  site-selective  Mott  transition)  and  the  demonstration  that  this  mechanism  ac¬ 
counts  for  the  main  features  of  the  previously  mysterious  metal-insulator  phase  dia¬ 
gram  of  the  rare-earth  nickelates. 

•  the  ability  to  design  new  Mott  insulators  via  heterostructuring. 

•  A  fully  charge  self-consistent  implementation  of  DFT-I-DMFT,  which  allows  us  to 
compute  total  energies.  This  major  technical  development  allowed  us  to  compute  the 
pressure  vs.  rare  earth  phase  diagram  in  the  nickelates. 
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1  Overview 

This  grant  funded  one  postdoctoral  researcher.  The  broad  goal  of  the  research  was  to  under¬ 
stand  the  degree  to  which  many-body  physics  could  be  manipulated  by  materials  design,  in 
particular  by  growth  of  materials  in  appropriately  chosen  atomic-thickness  heterostructure 
conhgurations.  Our  work  was  based  on  combining  the  most  advanced  many-body  tech¬ 
niques  available  to  treat  strongly  correlated  electron  systems  with  an  unprecedented  level 
of  physical  and  chemical  of  realism.  The  initial  focus  was  to  demonstrate  a  theoretical 
basis  for  control  of  one  fundamental  correlated  electron  property,  the  orbital  polarization 
(relative  occupancy  of  different  transition  metal  d  orbitals).  In  the  course  of  the  investiga¬ 
tion  we  uncovered  fundamental  difficulties  with  the  conventionally  accepted  theory  of  the 
correlation-driven  metal-insulator  transition,  which  we  investigated.  We  also  uncovered  a 
new  kind  of  correlation-driven  metal-insulator  transition  and  showed  that  our  mechanism 
accounts  for  the  previously  mysterious  properties  of  the  rare  earth  nickelates.  Addition¬ 
ally,  we  demonstrated  the  ability  to  design  new  Mott  insulators  via  heterostructured  double 
perovskites.  Finally,  we  succeeded  in  implementing  a  successful  fully  charge  self-consistent 
implementation  of  DFT-I-DMFT  which  allows  for  total  energy  calculations.  Our  work  has  led 
to  three  papers  fully  supported  by  the  grant  [1,  2,  3]  and  three  partially  supported[4,  22,  23]. 
Below  we  review  the  major  hndings. 

2  Accomplishments  in  previous  funding  period 

2.1  Orbital  Polarization:  the  case  of  Nickelate  heterostructures 

This  work[l,  2]  had  a  general  and  a  specihc  motivation.  The  general  motivation  is  to  es¬ 
tablish  the  body  of  theoretical  results  needed  for  the  broad  program  of  controlling  many 
body  physics  by  heterostructuring.  To  begin  work  on  this  large-scale  project  we  selected 
for  study  a  basic  property  (the  orbital  polarization  or  differential  occupancy  of  transition 
metal  d-orbitals)  which  is  fundamental  to  many-body  physics.  The  specihc  motivation  was 
a  proposal,  put  forth  by  Chaloupka  et  al[5]  and  supported  by  calculations  of  Hansmann 
et  al[6],  that  with  appropriate  heterostructuring  one  could  control  the  orbital  polarization 
to  the  extent  of  giving  a  nickelate  material  (two-sheeted  Fermi  surface)  a  cuprate-like  elec¬ 
tronic  structure,  thereby  producing  a  new  class  of  high  temperature  superconductors.  Our 
results  revised  the  conventional  picture  in  two  important  ways.  First,  we  showed  via  density 
functional  calculations  [1]  the  incorrectness  of  the  standard  belief  that  the  only  important 
effect  of  heterostructuring  was  geometric  conhnement.  In  fact,  the  chemical  composition  of 
the  supposed  electrically  inert  counterlayers  had  a  non-trivial  effect,  via  shifts  of  the  ligand 
energy  levels.  Second  we  showed  via  DFT-I-DMFT  calculations  [2]  that  for  ’late’  transition 
metal  oxides  such  as  the  lanthanum  nickelates  the  previously  accepted  theoretical  analysis 
was  irrelevant  because  it  neglected  the  physically  crucial  consequences  of  charge  transfer 
between  transition  metal  (here,  Ni)  and  ligand  (here,  O)  ions.  Taking  charge  transfer  and 
Hunds  coupling  physics  correctly  into  account  reveals  that  the  orbital  polarization  predicted 
by  Chaloupka  and  studied  (in  a  model  which  did  not  include  the  ligands)  by  Hansmann  et 
al  simply  does  not  occur  in  the  lanthanum  nickelate  materials;  indeed  when  the  correct 
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physics  is  incorporated  into  the  calculation  one  hnds  that  correlation  effects  act  to  decrease 
the  magnitude  of  the  orbital  polarization  (see  figure  1).  Very  recent  resonant  scattering 
experiments  appear  to  confirm  our  predictions  [7]. 


Figure  1:  The  Fermi  surface  of  nickelate  heterostructures  (layers  of  LaNiOa  alternating 
with  highly  insulating  counterlayers  along  the  perovskite  (111)  direction  using  DFT  and 
DFT+DMFT.  Panels  (a,b)  our  results  [2]  including  covalency  effects;  comparison  of  panels 
(a)  (DFT)  and  (b)  (DFT+DMFT)  shows  an  interaction-driven  decrease  in  orbital  polar¬ 
ization,  visible  here  as  an  increase  in  the  size  of  the  central  Fermi  surface.  Panels  (c,d): 
analogous  calculations [6]  but  with  many-body  physics  modelled  via  a  Hubbard  Hamiltonian 
which  does  not  include  covalency  effects.  Tendency  (in  this  approximation)  of  interactions 
to  increase  orbital  polarization  shown  by  correlation-driven  disappearance  of  central  pocket. 

Our  work  was  pivotal  both  in  clarifying  what  can  be  achieved  in  this  important  and 
widely  studied  system  and  more  broadly  in  revealing  importance  of  covalency  physics  pre¬ 
viously  neglected  in  many  studies  of  transition  metal  oxides. 

2.2  Covalency  and  the  metal-insulator  phase  diagram  in  transi¬ 
tion  metal  oxides 

The  work  described  in  the  previous  section  highlighted  the  importance  of  covalency  effects 
in  nickelates  and  raised  the  broader  question  of  the  role  of  covalency  effects  across  the  full 
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family  of  transition  metal  oxides.  The  large  and  growing  field  of  oxide  electronics  has  used 
the  Hubbard  model  (which  neglects  covalency  effects)  as  a  paradigm  for  strong  correlation 
physics,  and  it  is  important  to  know  how  covalency  effects  change  the  picture.  We  resolved 
these  issues[4],  and  we  describe  the  main  points  below. 

The  fundamental  parameter  controlling  covalency  is  the  energy  of  the  correlated  (e.g. 
transition  metal  d)  states  relative  to  the  other  electronic  states  (for  example  the  oxygen  p). 
Theoretical  analysis  raises  a  deep  question  in  materials  theory.  Methodologies  for  treating 
electronic  correlations  beyond  density  functional  theory  require  the  identihcation  of  a  set 
of  orbitals  whose  correlations  are  to  be  treated  more  accurately  (e.g.  the  d-like  orbitals 
centered  on  the  transition  metal  sites  of  transition  metal  oxides  or  the  /-like  states  centered 
on  the  actinides  in  heavy  Fermion  systems)  and  a  prescription  for  embedding  these  correlated 
orbitals  in  the  full  electronic  structure.  A  crucial  issue  in  the  embedding  prescription  is  the 
’double  counting  correction’  which  removes  the  correlation  contributions  from  the  single¬ 
particle  energies  obtained  in  density  functional  theory  and  thus  adjusts  the  energy  of  the 
correlated  orbitals  relative  to  other  electronic  states.  While  different  prescriptions  have  been 
proposed  [8,  9,  10],  there  is  no  broadly  agreed  upon  and  theoretically  justihed  expression 
for  the  double  counting  correction,  and  yet  because  it  controls  covalency  it  plays  a  crucial 
role  in  beyond  DFT  electronic  structure  methods.  Lack  of  theoretical  control  over  this  issue 
is  a  fundamental  limitation  of  modern  theories  of  correlated  electrons  in  solids. 


Figure  2:  Metal-insulator  phase  diagram  from  Ref.  [4]  indicating  limit  of  stability  of  metallic 
phase  Nd,c2  (solid  symbols),  limit  of  stability  of  insulating  phase  N^^ci  (open  symbols)  along 
with  vertical  lines  indicating  band  theory  estimates  of  Nd  for  model  of  high  Tc  cuprates 
(black  and  red  squares,  and  black  dashed  line)  and  nickelates  (blue  and  green  triangles,  and 
blue  dotted  line).  Note  that  in  single-site  dynamical  mean  held  theory  the  metal  insulator 
transition  is  hrst  order  and  is  characterized  by  two  spinodal  lines:  f/ci(A^ci)  where  the 
insulating  phase  loses  stability  and  a  slightly  different  value  Uc2{Nc2)  at  which  the  metallic 
phase  loses  stability  [11].  We  present  both  lines  here,  but  remark  that  the  Uc2  lines  have  a 
strong  temperature  dependence  and  are  difficult  to  compute  accurately. 


As  a  hrst  step  towards  an  improved  understanding  of  this  issue  we  have  investigated  [4] 
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different  ways  of  parameterizing  covalency  (and  therefore  the  double  counting  correction) 
and  of  comparing  the  results  to  experiment.  We  based  our  investigations  on  underlying  elec¬ 
tronic  structures  obtained  from  density  functional  band  theory,  with  the  relevant  correlated 
and  ligand  orbitals  dehned  by  a  maximally  localized  Wannier  function  construction,  and 
used  single-site  dynamical  mean  held  theory  (which  instantiates  the  canonically  accepted 
physics  of  the  Mott  metal-insulator  transition)  to  treat  the  correlation  effects.  We  focussed 
on  the  high-Tc  cuprates  (an  electronically  two  dimensional  system  with  one  relevant  cor¬ 
related  orbital)  and  the  rare  earth  nickelates  (electronically  three  dimensional,  with  two 
relevant  correlated  orbitals).  We  found  that  the  double  counting  and  covalency  effects  are 
most  usefully  parameterized  in  terms  of  the  occupancy  of  the  transition  metal  d-orbital. 
Different  theoretical  models  (involving  different  strengths  of  the  beyond  density  functional 
theory  interactions  and  different  double  counting  corrections)  which  yield  the  same  yield 
essentially  the  same  physics,  and  in  particular  the  metal-insulator  phase  diagram  takes  the 
nearly  universal  form  shown  in  the  right-hand  panel  of  Fig.  2.  (Note  that  while  the  precise 
value  of  of  course  depends  on  the  precise  dehnition  of  orbitals,  all  dehnitions  based  on 
well  localized  orbitals  give  very  similar  results).  A  particularly  striking  hnding  common 
to  both  systems  is  that  within  the  single-site  DMFT  approximation,  insulating  behavior  is 
found  only  for  a  very  narrow  range  of  d-occupancy,  corresponding  to  weak  covalency  effects, 
irrespective  of  the  Coulomb  repulsion.  Further,  while  as  remarked  above  the  double  count¬ 
ing  (and  hence  the  true  magnitude  of  covalency  effects)  is  not  known,  the  values  required 
to  make  the  theoretical  models  insulating  are  so  far  from  the  band  theory  values  or  litera¬ 
ture  estimates  from  various  experiments  that  it  is  highly  unlikely  that  within  the  single  site 
DMFT  approximation  and  therefore  the  canonical  Mott  picture  is  unable  to  account  for  the 
insulating  behavior  of  these  materials.  We  resolve  this  apparent  paradox  in  the  next  section 
where  we  describe  the  site-selective  Mott  transition. 

2.3  Metal  insulator  transition  in  rare  earth  nickelates. 

In  this  work,  we  resolved  the  nature  of  the  insulating  state  in  the  rare-earth  nickelates [3]. 
The  rare  earth  nickelates  have  emerged  as  one  of  the  key  paradigm  systems  for  oxide  elec¬ 
tronics  [12,  13,  14]  and  crucial  goal  of  current  research  is  to  understand  the  metal  insulator 
transitions  which  occur  as  a  function  of  chemical  composition,  him  thickness  and  imposed 
strain.  A  necessary  hrst  step  is  a  theory  for  the  metal-insulator  transition  occurring  in  the 
bulk  compounds.  The  investigations  outlined  above  show  that  this  transition  is  unlikely  to 
be  due  to  conventional  Mott  physics.  Experimental  work  [15]  has  showed  that  the  transition 
is  associated  with  a  structural  change  in  which  Ni-0  octahedra  become  inequivalent,  with 
one  having  a  decreased  and  the  other  an  increased  mean  Ni-0  bond  length.  These  results 
have  been  interpreted  as  evidence  for  a  crucial  role  of  charge  ordering  [16,  15,  17,  18,  19]  but 
charge  ordering  in  the  naive  sense  of  a  change  in  the  Ni  valence  between  the  two  sublattices 
is  inconsistent  with  the  large  repulsive  d-d  interaction  U  expected  on  the  Ni  site  and  appears 
to  be  ruled  out  by  very  recent  soft  X-ray  resonant  diffraction  data  [20].  Our  calculations 
also  reveal  negligible  charge  ordering.  Antiferromagnetic  order  can  drive  an  insulating  state 
and  induce  a  disproportionation  as  a  second-order  effect  [21];  while  this  mechanism  may  be 
relevant  to  NdNiOs  it  fails  to  account  for  the  paramagnetic  insulating  state  observed  in  Sm, 
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Eu,  Y,  and  Lu  compounds.  The  cause  of  the  transition  has  been  mysterious  but,  we  believe, 
our  new  results  solve  the  problem. 


b 


Figure  3:  Left  Panel  The  spectral  function  is  projected  onto  the  Wannier  basis  correspond¬ 
ing  to  the  Ni  eg  states  on  Nii  (dashed,  red  on-line)  and  Ni2  (solid,  blue  on-line)  and  displayed 
in  the  near  Fermi-surface  frequency  regime,  (a)  DFT  spectral  function,  (b)  DFT-I-DMFT 
spectral  function  obtained  at  U=5eV,  J=leV,  Nd=8.0,  and  T=116K.  Right  Panel:  The 
local  magnetic  susceptibility  x(i^  =  0)  in  a  paramagnetic  state  as  a  function  of  temperature 
is  obtained  using  the  DFT-I-DMFT  calculation.  Both  the  Nii  (large  Ni-0  bond  length,  red 
circle)  and  the  Ni2  (small  bond  length,  blue  square)  are  compared.  In  the  inset,  the  static 
magnetic  moments  as  a  function  of  temperature  in  Nii(red  circle)  and  Ni2(blue 

square)  are  computed  within  DFT-I-DMFT  in  a  ferromagnetic  state. 

We  studied  the  end  member  of  the  rare-earth  nickelate  series,  LuNiOs,  using  state  of  the 
art  density  functional  theory  plus  dynamical  mean  field  theory  (DFT-I-DMFT)  methods  to 
investigate  the  interplay  between  correlation  physics  and  structural  distortions.  Our  cal¬ 
culations  used  the  real  crystal  structure,  including  the  two  sublattice  difference  in  Ni  —  O 
octahedral  size.  We  explored  various  values  of  the  double  counting  correction  but  found  that 
the  key  results  do  not  depend  on  the  double  counting  correction  in  any  signihcant  way.  Our 
main  results  are  shown  in  Fig.  3.  The  left  panel  compares  the  d-electron  spectral  function 
for  LuNiOs  obtained  from  a  density  functional  calculation  (Fig.  3a)  to  that  obtained  from 
our  DFT-I-DMFT  calculation  (Fig.  3b).  We  see  that  while  density  functional  theory  pre¬ 
dicts  that  the  system  is  metallic  (no  gap  at  the  chemical  potential)  despite  the  bond  length 
disproportionation,  the  correlation  effects  captured  by  DMFT  drive  the  system  into  an  insu¬ 
lating  state  with  a  gap  of  ~  200meV  (the  precise  gap  value  depends  on  the  double  counting 
prescription).  Further,  although  as  many  workers  have  noted  the  disproportionation  im¬ 
plies  a  difference  in  charge  between  the  two  Ni  sites,  we  find  that  the  total  d  occupancy 
is  8.22  ±  0.08  for  the  DFT  spectra  and  8.0  ±  0.06  for  the  DFT-I-DMFT  spectra  shown  in 
Fig.  3;  our  GGA-I-U  calculations  (not  shown)  give  an  even  smaller  difference  between  the 
two  Ni  sites.  We  therefore  conclude  that  the  charge  order,  while  of  course  nonvanishing  by 
symmetry,  is  not  physically  important  and  in  particular  is  not  the  cause  of  the  insulating 
behavior,  which  is  instead  due  to  a  novel  kind  of  insulating  state  which  we  have  uncovered. 
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To  demonstrate  the  new  features  of  the  insulating  state  found  in  our  DFT+DMFT  cal¬ 
culations  we  show  in  Fig.  3  the  temperature  dependence  of  the  local  magnetic  susceptibility. 
A  dramatic  difference  is  visible  between  the  two  sites.  The  local  susceptibility  of  the  Nii 
(long  bond)  site  has  a  well  dehned  1/T  behavior  indicating  well  dehned  and  long-lived  but 
thermally  fluctuating  magnetic  moments,  as  expected  in  a  paramagnetic  Mott  insulator  (in 
the  single-site  DMFT  approximation).  By  contrast,  the  local  susceptibility  of  the  Ni2  (short 
bond)  site  has  a  negligible  temperature  dependence,  indicating  a  lack  of  a  local  moment  on 
this  site.  For  this  reason,  we  identify  the  origin  of  the  insulating  state  as  a  site-selective 
Mott  transition.  The  Ni-0  covalency  plays  a  crucial  role  in  this  insulating  state.  The  Ni 
valence  is  near  d®  (spin  1  on  each  Ni),  and  there  is  correspondingly  one  hole  per  every  three 
oxygen.  The  Nii  site  displays  the  physics  of  a  S=1  Mott  insulator,  whereas  the  two  holes  on 
the  octahedron  surrounding  the  Ni2  site  bind  with  the  S'  =  1  on  the  Ni2  to  form  a  singlet; 
this  binding  energy  gives  the  excitation  gap  of  the  insulating  state. 


3  Progress  during  final  year 

Here  we  summarize  the  results  that  were  obtained  in  the  hnal  year  of  funding.  A  total  of 
two  publications  were  produced  during  this  time  period[22,  23]. 

3.1  Designing  a  S  =  1  Mott  insulator  via  heterostructuring 

The  results  of  this  section  were  published  in  Physical  Review  Letters [22].  In  this  work,  the 
goal  was  to  design  a  new  Mott  insulator  via  heterostructuring.  We  took  two  base  materials, 
LaNiOs  and  LaTiOs,  and  combined  them  in  a  heterostructure.  In  the  bulk,  LaNiOs  is  a 
metal  and  LaTiOa  is  a  Mott  insulator.  This  is  somewhat  easy  to  anticipate  given  that  Ti  is 
in  a  nominal  conhguration  and  Ni  is  in  a  nominal  dJ  conhguration.  Due  to  the  different 
relative  energies  of  the  valence  states  in  these  materials,  charge  will  nominally  flow  when 
they  are  intertwined  in  a  heterostructure.  As  one  might  expect,  the  Ti  d  states  will  be 
higher  than  the  Ni  d  states,  and  this  is  illustrated  schematically  in  hgure  4  (left  panel). 
Therefore,  it  is  expected  that  charge  may  flow  from  Ti  to  Ni,  converting  LaTiOs  to  a  band 
insulator  and  driving  Ni  to  a  2-\-  state.  Ni^+  may  be  expected  to  be  a  Mott  insulator,  as  is 
NiO.  These  simple  ideas  can  be  tested  in  actual  calculations  (see  hgure  4  right  panel).  As 
expected,  lining  up  the  upper  edge  of  the  oxygen  valence  band  indicates  that  the  Ti  will 
likely  donate  charge.  Indeed,  the  heterostructure  does  end  up  being  a  A  =  1  Mott  insulator. 
The  similarity  to  NiO  can  clearly  be  seen  in  the  density  of  states. 

This  work  demonstrates  that  simple  intuitive  reasoning  combined  with  complex  calcu¬ 
lations  can  straightforwardly  be  used  to  design  novel  Mott  insulators.  More  generally,  this 
system  can  be  viewed  as  a  double  perovskite,  which  refers  to  the  fact  that  we  have  two 
transition  metals.  We  believe  that  increasing  the  complexity  of  heterostructures  by  moving 
to  double  perovskites  will  open  up  many  new  opportunities. 


Project  Narrative:  Page  6 


C.A.  Marianetti  and  A.J.  Millis  -  Columbia  University 


Figure  4:  (Left  Panel  a)  Schematic  band  structures  of  com-  ponent  materials  LaTiOs  and 
LaNiOs-  The  dashed  purple  lines  are  the  Fermi  levels  for  the  two  materials.  LaTiOs  shows 
insulating  behavior  with  a  small  excitation  gap  set  by  Ti  d-d  transitions  and  a  wide  energy 
separation  between  Ti  d  states  and  O  p  states.  LaNiOS  exhibits  metallic  behavior  with 
strong  mixing  between  Ni  d  states  and  O  p  states.  The  red  arrow  highlights  the  direction 
of  charge  transfer  in  the  superlattice.  (Left  panel  b  )  Schematic  band  structure  of  the 
(LaTi03)i/(LaNi03)i  superlattice.  Ti  d  states  are  above  the  Fermi  level  (dashed  purple 
line).  Correlation  effects  split  Ni  d  states  into  lower  and  upper  Hubbard  bands,  separated 
by  Uatj  .  (Right  Panel)  Densities  of  states  for  majority  (above  axis)  and  minority  (below 
axis)  spins  of  super-  lattice  (upper  left)  and  reference  materials  NiO  (lower  left),  LaTiOs 
(upper  right;  zero  of  energy  is  shifted  so  that  oxygen  bands  align  with  those  of  LaNiOs  ) 
and  LaNiOs  (lower  right).  The  densities  of  states  are  obtained  using  DFT-I-U  calculations 
with  Utvi  =  6eld  and  U-ri  =  4eH. 

3.2  Computing  the  Pressure  vs.  Rare-Earth  Phase  digram  in  the 
Nickelates 

The  results  of  this  section  were  submitted  to  Physical  Review  Letters[23].  We  have  im¬ 
plement  a  fully  charge-self-consistent  density  functional  theory  plus  dynamical  mean  field 
theory  (DFT-I-DMFT)  methodology  which  enables  relatively  efficient  calculation  of  the  total 
energy  of  realistic  correlated  electron  systems.  The  DFT  portion  of  the  calculation  uses  a 
plane  wave  basis  set  within  the  project  augmented  wave  method;  enabling  study  of  systems 
with  large  unit  cells.  The  DMFT  part  uses  maximally  localized  Wannier  function  (MLWF) 
methods  to  define  the  correlated  orbitals,  enabling  a  computationally  efficient  interface  with 
the  DFT  calculation.  The  sensitivity  of  the  results  to  the  choice  of  double  counting  cor¬ 
rection  is  demonstrated  and  a  physically  reasonable  construction  is  motivated  and  verified 
by  comparison  to  experiment.  The  formalism  is  validated  by  application  to  the  rare-earth 
nickelate  family  of  materials.  The  metal-insulator  and  structural  phase  boundaries  are  de¬ 
termined  as  a  function  of  applied  pressure  and  choice  of  A-site  ion  and  are  found  to  be 
consistent  with  experiment.  The  DFT-|-f/  method  is  implemented  within  the  same  basis  set 
and  is  shown  to  grossly  overestimates  the  tendency  to  bond-disproportionation. 
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Our  total  energy  calculation  on  LuNiOa  is  illustrated  in  the  left  panel  of  Figure  5.  First, 
we  recall  that  DFT  completely  fails  to  to  predict  the  distortion  LuNiOs.  Alternatively, 
DFT+U  does  predict  a  transition  LuNiOs,  though  it  overpredicts,  while  it  also  erroneously 
predicts  a  transition  in  LaNiOs.  Our  DFT+DMFT  calculations  are  shown  to  successfully 
predict  the  distortion  in  LuNiOa,  though  the  magnitude  is  slightly  underpredicted.  Also 
shown  in  this  plot  is  the  effect  of  including  full  charge  self-consistency.  Somewhat  unexpect¬ 
edly,  the  effect  is  not  terribly  signihcant.  This  is  encouraging  in  the  sense  that  it  suggests 
that  non-charge-self-consistent  calculations  could  be  viable  for  total  energies. 

We  computed  the  pressure  phase  diagram  for  the  entire  rare-earth  nickelate  series,  and 
this  is  shown  in  the  right  panel  of  Figure  5.  The  black  line  corresponds  to  experimental 
measurements.  The  dotted  lines  correspond  to  standard  DFT-I-U  calculations,  which  are 
clearly  in  poor  agreement  with  experiment.  The  solid  lines  are  the  result  of  DFT-I-DMFT 
calculations,  and  they  are  in  excellent  agreement  with  experiment.  This  powerful  result 
demonstrates  that  DFT-I-DMFT  is  a  robust  method  for  computing  total  energies  in  strongly 
correlated  electron  systems. 


Figure  5:  (Left  Panel)  Total  energy  as  a  function  of  bond-length  difference  6a  for  LuNiOs 
calculated  using  fully  charge  self  consistent  DFT-I-DMFT  with  original  (squares,  black) 
and  modihed  (diamonds,  blue)  double  counting  correction  and  compared  to  DFT-I-DMFT 
energies  computed  using  the  DFT  charge  density  (N.C.S.C,  pentagons,  green)  and  to  en¬ 
ergies  obtained  from  the  DFT-|-t/  method  (circles,  red).  The  experimentally  determined 
values  {Sa  =  O.OSSA  (LuNiOs)  and  5a  =  0  (LaNiOs))  are  indicated  by  vertical  dashed  lines. 
The  parameters  are  f/=5.0eV,  L=1.0eV,  and  T=116K.  (Right  Panel)  Metal-insulator  and 
structural  phase  diagram  computed  using  DFT-I-DMFT  (solid  symbols  and  solid  lines)  as 
a  function  of  unit  cell  volume  for  the  series  of  rare  earth  ions  and  compared  to  results  of 
DFT-|-[/  calculations  (open  symbols,  dashed  lines)  and  to  experimental  data  (pentagons 
and  dashed  lines)  obtained  for  {Y ,Eu,Nd,Pr)'Ni03.  Vq  is  determined  as  the  calculated  am¬ 
bient  pressure  equilibrium  volume  for  each  material  while  the  tolerance  factor  is  determined 
from  the  distances  du-o  and  dm-o  as  dn-o/dNi-oV^-  The  parameters  for  the  DMFT 
calculations  are  r=116K,  U=5eV,  17'=4.8eV  and  L=leV. 
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